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Abstract 

The search for MSSM Higgs bosons will be an important goal at the LHC. We 
analyze the search reach of the CMS experiment for the heavy neutral MSSM Higgs 
bosons with an integrated luminosity of 30 or 60 fb^^. This is done by combining the 
latest results for the CMS experimental sensitivities based on full simulation studies 
with state-of-the-art theoretical predictions of MSSM Higgs-boson properties. The 
results are interpreted in MSSM benchmark scenarios in terms of the parameters tan [5 
and the Higgs-boson mass scale. Ma- We study the dependence of the 5 a discovery 
contours in the M^-tan (5 plane on variations of the other supersymmetric parameters. 
The largest effects arise from a change in the higgsino mass parameter /i, which enters 
both via higher-order radiative corrections and via the kinematics of Higgs decays into 
supersymmetric particles. While the variation of ^ can shift the prospective discovery 
reach (and correspondingly the "LHC wedge" region) by about Atan/3 = 10, we 
find that the discovery reach is rather stable with respect to the impact of other 
supersymmetric parameters. Within the discovery region we analyze the accuracy 
with which the masses of the heavy neutral Higgs bosons can be determined. We 
find that an accuracy of 1-4% should be achievable, which could make it possible in 
favourable regions of the MSSM parameter space to experimentally resolve the signals 
of the two heavy MSSM Higgs bosons at the LHC. 
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1 Introduction 



Identifying the mechanism of electroweak symmetry breaking will be one of the main goals of 
the LHC Many possibilities have been studied in the literature, of which the most popular 
ones are the Higgs mechanism within the Standard Model (SM) and within the Minimal 
Supersymmetric Standard Model (MSSM) [1]. Contrary to the case of the SM, in the MSSM 
two Higgs doublets are required. This results in five physical Higgs bosons instead of the 
single Higgs boson of the SM. These are the light and heavy CP-even Higgs bosons, h and 
H, the CP-odd Higgs boson. A, and the charged Higgs boson. The Higgs sector 

of the MSSM can be specified at lowest order in terms of the gauge couplings, the ratio 
of the two Higgs vacuum expectation values, tan/3 = V2/V1, and the mass of the CP-odd 
Higgs boson. Ma- Consequently, the masses of the CP-even neutral Higgs bosons and the 
charged Higgs boson are dependent quantities that can be predicted in terms of the Higgs- 
sector parameters. Higgs-phenomenology in the MSSM is strongly affected by higher-order 
corrections, in particular from the sector of the third generation quarks and squarks, so that 
the dependencies on various other MSSM parameters can be important. 

After the termination of LEP in the year 2000 (the final LEP results can be found in 
Refs. [2,3]), and the (ongoing) Higgs boson search at the Tevatron [4-6], the search will be 
continued at the LHC [7-9] (see also Refs. [10, 11] for recent reviews). The current exclusion 
bounds within the MSSM [3-5] and the prospective sensitivities at the LHC are usually dis- 
played in terms of the parameters Ma and tan/3 that characterize the MSSM Higgs sector 
at lowest order. The other MSSM parameters are conventionally fixed according to certain 
benchmark scenarios [12-14] . The most prominent one is the "m™^'' scenario" , which in the 
search for the light CP-even Higgs boson allows to obtain conservative bounds on tan /3 for 
fixed values of the top-quark mass and the scale of the supersymmetric particles [15]. Besides 
the "no-mixing scenario" , which is similar to the m™*^^ scenario, but assumes vanishing mix- 
ing in the stop sector, other CP-conserving scenarios that have been studied in LHC analyses 
(see e.g. Ref. [11]) are the "gluophobic Higgs scenario" and the "small aes" scenario [13]. 

For the interpretation of the exclusion bounds and prospective discovery contours in 
the benchmark scenarios it is important to assess how sensitively the results depend on 
those parameters that have been fixed according to the benchmark prescriptions. While in 
the decoupling limit, which is the region of MSSM parameter space with Ma ^ Mz, the 
couplings of the light CP-even Higgs boson approach those of a SM Higgs boson with the 
same mass, the couplings of the heavy Higgs bosons of the MSSM can be sizably affected 
by higher-order contributions even for large values of Ma- The kinematics of the heavy 
Higgs-boson production processes, on the other hand, is governed by the parameter Ma, 
since in the region of large Ma the heavy MSSM Higgs bosons are nearly mass-degenerate. 
Ma ~ Mh ~ Mh±. In Ref. [14] it has been shown that higher-order contributions to the 
relation between the bottom-quark mass and the bottom- Yukawa coupling have a dramatic 
effect on the exclusion bounds in the My^-tan (3 plane obtained from the hhcj), cf) ^ hh channel 
at the Tevatron. 

In this article we investigate how the 5 a discovery regions in the M^-tan (3 plane for the 
heavy neutral MSSM Higgs bosons (a corresponding analysis for the charged Higgs-boson 

^We focus in this paper on the case without explicit CT'-violation in the soft supersymmetry-breaking 
terms. 
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search will be presented elsewhere) obtainable with the CMS experiment at the LHC depend 
on the other MSSM parameters. For the experimental sensitivities achievable with CMS we 
use up-to-date results based on full simulation studies for 30 or 60 fb~^ (depending on the 
channel) [9]. This information is combined with precise theory predictions for the Higgs- 
boson masses and the involved production and decay processes incorporating higher-order 
corrections at the one-loop and two-loop level. In our analysis we investigate the impact on 
the discovery reach arising both from higher-order corrections and from possible decays of 
the heavy Higgs bosons into supersymmetric particles^] 

The search for the heavy neutral MSSM Higgs bosons at the LHC will mainly be pursued 
in the b quark associated production with a subsequent decay to r leptons [7-9]. In the region 
of large tan /3 this production process benefits from an enhancement factor of tan^ f3 compared 
to the SM case. The main search channels ar^ (here and in the following (p denotes the two 
heavy neutral MSSM Higgs bosons, (p = H,A): 



6^ r+r" ^2jets (1) 

bb(f), r+r" ^ + jet (2) 

b T^T' ^ e + jet (3) 

b r+r" ^ e + /i . (4) 



For our numerical analysis we use the program FeynHiggs [19-22]. We study in particular 
the dependence of the "LHC wedge" region, i.e. the region in which only the light CP-even 
MSSM Higgs boson can be detected at the LHC at the 5 a level, on the variation of the 
higgsino mass parameter /i. The dependence on n enters in two different ways, on the one 
hand via higher-order corrections affecting the relation between the bottom mass and the 
bottom Yukawa coupling, and on the other hand via the kinematics of Higgs decays into 
supersymmetric particles. We analyze both effects separately and discuss the possible impact 
of other supersymmetric parameters. 

Our results for the discovery reach of the heavy neutral MSSM Higgs bosons extend the 
known results in the literature in various ways. In comparison with Refs. [23,24], where the 
prospective 5a discovery contours for CMS in the M^-tan/? plane of the m^^^ benchmark 
scenario were given for three different values of /x, the results in the present paper are based 
on full simulation studies and make use of the most up-to-date CMS tools for triggering 
and event reconstruction. Furthermore, in the analysis of Refs. [23,24] relevant higher- 
order corrections, in particular those depending on (see Sect. 12.21 below), have been 
neglected. The effects induced by the Af, corrections have been investigated in Ref. [14], 
where the results were obtained by a simple rescaling of the experimental results given in 
Refs. [7, 23-25] . Our present analysis, on the other hand, makes use of the latest CMS studies 
and provides a separate treatment of the different r final states, channels 

As a second step of our analysis we investigate the experimental precision that can be 
achieved for the determination of the heavy Higgs-boson masses in the discovery channels ([1])- 



^We restrict our analysis to the impact of supersymmetric contributions. For a discussion of uncertainties 
related to parton distribution functions, see e.g. Ref. [16]. 

■^In our analysis we do not consider diffractive Higgs production, pp p ® H ® p [17]. For a detailed 
discussion of the search reach for the heavy neutral MSSM Higgs bosons in diffractive Higgs production we 
refer to Ref. [18]. 
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(jl]). We discuss the prospective accuracy of the mass measurement in view of the possibihty 
to experimentally resolve the signals of the heavy neutral MSSM Higgs bosons. 

The paper is organized as follows: Sect. [2] introduces our notation and gives a brief sum- 
mary of the most relevant supersymmetric radiative corrections to the Higgs-boson masses, 
production cross sections and decay widths at the LHC The relevant benchmark scenarios 
are briefly reviewed. In Sect. [3] the experimental analysis is described. The results for the 
variation of the 5 a discovery contours, obtainable at CMS with 30 or 60 fb~^ are given 
in Sect. HI where we also discuss the achievable experimental precision in the Higgs mass 
determination. The conclusions can be found in Sect. [H 



2 Phenomenology of the MSSM Higgs sector 
2.1 Notation 

The MSSM Higgs sector at lowest order is described in terms of two independent parameters 
(besides the SM gauge couplings): tan/5 = V2/V1, the ratio of the two vacuum expectation 
values, and M4, the mass of the CP-odd Higgs boson A. Beyond the tree- level, large 
radiative corrections can occur from the t/i sector, and for large values of tan/3 also from 
the b/b sector. 

Our notations for the scalar top and scalar bottom sector of the MSSM are as follows: 
the mass matrices in the basis of the current eigenstates t^, t/j and &l, 6_r are given by 

A.2 _ / M?+m? + COS 2/? (i-|4)M| m,X, \ 

~ 1, mtXt M^ + m2 + |cos2/3 4M| j ' 

A.2 _ /M?+m2 + cos2/?(-i + |4)M| m,X, \ 



where 



mtXt = mt{At — fi cot P ) , rubXh = nib {A^ — fi tan P). (7) 



Here Mq, M^^ and M^^ are the diagonal soft SUSY-breaking parameters. At denotes the 
trilinear Higgs-stop coupling, A^ denotes the Higgs-sbottom coupling, and fi is the higgsino 
mass parameter. 

For the numerical evaluation, it is often convenient to choose 



Mq = il^^ = M^^ =: M 



SUSY- 



Concerning analyses for the case where M^^ ^ Mq 7^ ^6^5 see e.g. Refs. [20,26,27]. It 
has been shown that the upper bound on the mass of the light CP-even Higgs boson, M/j, 
obtained using eq. ([8]) is the same as for the more general case, provided that Msusy is 
identified with the heaviest mass of Mq, M^^, M^^ [20]. 

Accordingly, the most important parameters entering the Higgs-sector predictions via 
higher-order corrections are rrit, MgusY, Xt, Xj, and fi (see also the discussion in Sect. 12.2.21 
below). The Higgs-sector observables furthermore depend on the SU(2) gaugino mass param- 
eter, M2, the U(l) parameter Mi and the gluino mass, rrig (the latter enters the predictions 
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for the Higgs-boson masses only from two-loop order on). In numerical analyses the U(l) 
gaugino mass parameter, Mi, is often fixed via the GUT relation 

Ml = 14^2. (9) 

We will briefly comment below on the possible impact of complex phases entering the Higgs- 
sector predictions via higher-order contributions. 

2.2 Higher-order corrections in the Higgs sector 

In the following we briefly summarize the most important higher-order corrections affecting 
the observables in the MSSM Higgs-boson sector. As mentioned above, we focus on the 
MSSM with real parameters. For our numerical analysis we use the program FeynHiggs [19- 
220, which incorporates a comprehensive set of higher-order results obtained in the Feynman- 
diagrammatic approach [20-22,28-30]. 

2.2.1 Higgs-boson propagator corrections 

Higher-order corrections to the Higgs-boson masses and the wave function normalization 
factors of processes with external Higgs bosons arise from Higgs-boson propagator-type con- 
tributions. These corrections furthermore contribute in a universal way to all Higgs-boson 
couplings. For the propagator-type corrections in the MSSM the complete one-loop re- 
sults [31-34], the bulk of the two-loop contributions [20,27-29,35-39] and even leading 
three-loop corrections [40] are known. The remaining theoretical uncertainty on the light 
CP-even Higgs-boson mass has been estimated to be below ~ 3 GeV [21,41]. The by far 
dominant contribution is the 0{at) term due to top and stop loops {at = /i^/(47r), where ht 
denotes the top-quark Yukawa coupling). Effects of 0{ab) can be important for large values 
of tan p. 

2.2.2 Corrections to the relation between the bottom-quark mass and the bot- 
tom Yukawa coupling 

Concerning the corrections from the bottom/sbottom sector, large higher-order effects can 
in particular occur in the relation between the bottom-quark mass and the bottom Yukawa 
coupling (which controls the interaction between the Higgs bosons and bottom quarks as 
well as between the Higgs and scalar bottoms), hb, for large values of tan/5. At lowest order 
the relation reads rrib = hbVi. Beyond the tree level large radiative corrections proportional 
to hbV2 are induced, giving rise to tan /5-enhanced contributions [36-38,42]. At the one-loop 
level the leading terms proportional to V2 are generated either by gluino-sbottom one-loop 
diagrams of 0{as) or by chargino-stop loops of 0{at). 

The leading one-loop contribution in the limit of Msusy 3> rrit and tan /? ^ 1 takes 
the simple form [36] 

Ab = nig fi tan P x /(mi , mt , m^) + — At tan/3 x I{mi ,mi^) , (10) 
The code can be obtained from www.feynhiggs.de . 
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where the function I is given by 

na,b,c) = _ ^.)(^. ' _ [-'b' log ^ + b^c^ log ^ + c^a^ log ^) (11) 

1 



The leading contribution can be resummed to all orders in the perturbative expansion [BO- 
SS]. This leads in particular to the replacement 

where nib denotes the running bottom quark mass including SM QCD corrections. For the 
numerical evaluations in this paper we choose mj, = rnb{mt) ~ 2.97 GeV. 

The Afe corrections are numerically sizable for large tan /3 in combination with large values 
of the ratios of /im^/MgusY /^^t/^susY- Negative values of lead to an enhancement of 
the bottom Yukawa coupling as a consequence of eq. fll2l) (for extreme values of /i and tan (3 
the bottom Yukawa coupling can even acquire non-perturbative values when Af, — — 1), 
while positive values of Af, give rise to a suppression of the Yukawa coupling. Since a change 
in the sign of n reverses the sign of Af,, the bottom Yukawa coupling can exhibit a very 
pronounced dependence on the parameter /i. 

For large values of tan/5 the correction to the production cross sections of the Higgs 
bosons H and A induced by Af, enters approximately like tan^/?/(l + A;,)^, giving rise 
to potentially large numerical effects. In the case of the subsequent Higgs-boson decay 
r+r^, however, the A;, corrections in the production and the decay process cancel 
each other to a large extent. The residual A^ dependence of a{bb(f)) x BR(0 r+r^) is 
approximately given by tan^ /5/((l + Af,)^ + 9), which has a much weaker A^ dependence (see 
Ref. [14] for a more detailed discussion). 

In the numerical analysis below the A^ corrections, which have been discussed in this 
section in terms of simple approximation formulae, will be supplemented by other higher- 
order corrections as implemented in the program FeynHiggs (and possible decay modes into 
supersymmetric particles are taken into account). Higher-order corrections to Higgs decays 
into r+r" within the SM and MSSM have been evaluated in Refs. [34,43]. 



2.2.3 Corrections to the Higgs production cross sections 

For the prediction of Higgs-boson production processes at hadron colliders SM-type QCD 
corrections in general play an important role. The SM predictions for the process bb (p+X 
at the LHC are far advanced. In the five-flavor scheme the SM cross section is known at 
NNLO in QCD [44]. The cross section in the four-flavor scheme is known at NLO [45,46]. 
Results obtained in the two schemes have been shown to be consistent [47-49] (see also 
Refs. [48,50] and Refs. [45,46] for results with one and two final-state 6-quarks at high-pj-, 
respectively). 

The predictions for the bb ^ (p + X cross sections in the MSSM have been obtained 
with FeynHiggs [19-22]. The FeynHiggs implementatiorl^l is based on the state-of-the-art 

■''The inclusion of the charged Higgs production cross sections is planned for the near future. 
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SM prediction, namely the NNLO result in the five-flavor scheme [44] using MRST2002 
parton distributions at NNLO [51], with the renormalization scale set equal to Mjjsm and 
the factorization scale set equal to Mhsm/4. In order to obtain the MSSM prediction the 
SM cross section is rescaled with the ratio of the partial widths in the MSSM and the SM, 

r(0 bb)sM 

The evaluation of the partial widths incorporates one-loop SM QCD and SUSY QCD correc- 
tions, as well as (in the SUSY case) the resummation of all terms of 0((astan/3)"') [34,37,43] 
and the proper normalization of the external Higgs bosons as discussed in Refs. [22, 52]. Since 
the approximation of rescaling the SM cross section with the ratio of partial widths does not 
take into account the MSSM-specific dynamics of the production processes, the theoretical 
uncertainty in the predictions for the cross sections will in general be somewhat larger than 
for the decay widths. It should be noted that in comparison with other approaches for treat- 
ing the SM and SUSY contributions, for instance the program HQQ [53], sizable deviations 
can occur as a consequence of differences in the scale choices and the inclusion of higher-order 
corrections. 



2.3 The m™^^ and no-mixing benchmark scenarios 

While the phenomenology of the production and decay processes of the heavy neutral MSSM 
Higgs bosons at the LHC is mainly characterised by the parameters M4 and tan (3 that govern 
the Higgs sector at lowest order, other MSSM parameters enter via higher-order contribu- 
tions, as discussed above, and via the kinematics of Higgs-boson decays into supersymmetric 
particles. The other MSSM parameters are usually fixed in terms of benchmark scenarios. 
The most commonly used scenarios are the "m™^^" and "no-mixing" benchmark scenar- 
ios [12-14]. According to the definition of Ref. [13] the m™'^^ scenario is given by 

mr'' : MsusY = 1000 GeV, Xt = 2MsusY, A = A, 

/i = 200 GeV, M2 = 200 GeV, rrig = 0.8 MgusY • (14) 

The no-mixing scenario differs from the m™'^^ scenario only in that it has vanishing mixing 
in the stop sector and a larger value of Msusy 

no- mixing: Msusy = 2000 GeV, Xt = 0, Ab = At, 

fi = 200 GeV, M2 = 200 GeV, m-g = 0.8 Msusy • (15) 

The value of the top-quark mass in Ref. [13] was chosen according to the experimental central 
value at that time. For our numerical analysis below, we use the value, mt = 171.4 GeV [54j^. 

In Ref. [14] it was suggested that in the search for heavy MSSM Higgs bosons the m^^^ 
and no-mixing scenarios, which originally were mainly designed for the search for the light 
CP-even Higgs boson h, should be extended by several discrete values of /i, 

= ±200, ±500, ±1000 GeV . (16) 

^ Most recently the central experimental value has shifted to mt = 170.9 ± 1.8 GeV [55]. This shift has 
a negligible impact on our analysis. 
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As discussed above, the variation of fi in particular has an impact on the correction A;,, 
modifying in this way the bottom Yukawa couphng. For very large values of tan/3 and 
large negative values of fi the bottom Yukawa coupling can be so much enhanced that a 
perturbative treatment is no longer possible. We have checked that in our analysis of the 
LHC discovery contours the bottom Yukawa coupling stays in the perturbative regime, so 
that all values of /i down to /i = —1000 GeV can safely be inserted. 

The variation of the parameter /i also modifies the mass spectrum and the couplings in 
the chargino and neutralino sector of the MSSM. Besides the small higher-order corrections 
induced by loop diagrams involving charginos and neutralinos, a change in the mass spectrum 
of the chargino and neutralino sector can have an important effect on Higgs phenomenology 
because decay modes of the heavy neutral MSSM Higgs bosons into charginos and neutralinos 
open up if the supersymmetric particles are sufficiently light (the mass spectrum in the m™'^^ 
and no-mixing scenarios respects the limits from direct searches for charginos at LEP [56] 
for all values of /i specified in eq. f|T6|) ). 

Differences between the m™^^ and no-mixing scenarios in the searches for heavy neutral 
MSSM Higgs bosons are induced in particular by a difference in the A;, correction. While in 
the m™^"" scenario both the 0{as) and 0{at) contributions to A^ can be sizable, see eq. (fTOj) . 
in the no-mixing scenario the 0{at) contribution is very small because At is close to zero in 
this case. The larger value of Msusy in the no-mixing scenario gives rise to an additional 
suppression of |Ah| compared to the m™^"^ scenario. 

3 Experimental analysis 

In this section we briefly review the recent CMS analysis of the (p — >■ r+r^ channel, see 
Ref. [9], yielding the number of events needed for a 5 a discovery (depending on the mass 
of the Higgs boson). The analysis was performed with full CMS detector simulation and 
reconstruction for the following four flnal states of di-r-lepton decays: r+r^ jets [57], 
r+r" ^ e + jet [58], r+r" ^ /i jet [59] and r+r" ^ e + jj [60]. 

The Higgs-boson production in association with b quarks, pp — > bb(f), has been selected 
using single 6-jet tagging in the experimental analysis. The kinematics of the gg bb(j) 
production process (2 3) was generated with PYTHIA [61]. It has been shown that in 
this way the NLO kinematics is better reproduced than using the PYTHIA gb b(f) process 
(2 — > 2) [62]. The backgrounds considered in the analysis were QCD muli-jet events (for the 
rr — > jets mode), ti, bb, Drell-Yan production of Z, 7*, VT+jet, Wt and rrbb. All background 
processes were generated using PYTHIA, except for T~^T~bb, which was generated using 
CompHEP [63]. 

The results for the various channels, eqs. ([I]) - dl]), are given in Tabs. [T]-|H For every 
Higgs-boson mass point studied we show the number of signal events needed for 5 cr discovery, 
Ns, the total experimental selection efficiency, £exp, and the ratio of the di-r mass resolution 
to the Higgs-boson mass, Rm^- The last row in Tabs. [T]-|l] shows the expected precision of 
the Higgs-boson mass measurement, evaluated as explained below, for parameter points on 
the 5 a discovery contour. Detector effects, experimental systematics and uncertainties of 
the background determination were taken into account in the evaluation of the A^^. These 
effects reduce the discovery region in the M^-tan j3 plane as shown in previous analyses [9] 
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(j) T+T- jets, 60 fb"^ 


Ma [GeV] 


200 


500 


800 


Ns 


63 


35 


17 


^cxp 


2.5 X 10-^ 


2.4 X 10-3 


3.6 X 10-3 




0.176 


0.171 


0.187 


AM^/M^ [%] 


2.2 


2.8 


4.5 



Table 1: Required number of signal events, Ns, with £ = 60 fb^ for a 5o" discovery in 
the channel (p ~^ t^t~ — > jets. Furthermore given are the total experimental selection 
efficiency, ^exp, the ratio of the di-r mass resolution to the Higgs-boson mass, Rm^, and the 
expected precision of the Higgs-boson mass measurement, AM^/M^, obtainable from Ns 
signal events. 



T+T- ^ e + jet, 30 fb-^ 


Ma [GeV] 


200 


300 


500 


Ns 


72.9 


45.5 


32.8 


^exp 


3.0 X 10-3 


6.4 X 10-3 


1.0 X 10-2 


Rm^ 


0.216 


0.214 


0.230 


AM^/M^ [%] 


2.5 


3.2 


4.0 



Table 2: Required number of signal events, Ns, with £ = 30 fb for a 5 cr discovery in the 
channel —>■ t^t^ e+ jet. The other quantities are defined as in Tab. [TJ 



r+r- ^ ^1+ jet, 30 fb-^ 


Ma [GeV] 


200 


500 


Ns 


79 


57 


^exp 


7.0 X 10-3 


2.0 X 10-2 


Rm^ 


0.210 


0.200 


AM^/M^ [%] 


2.4 


2.6 



Table 3: Required number of signal events, Ns, with £ = 30 fb for a 5 cr discovery in the 
channel —>■ t^t^ —>■ fi + jet. The other quantities are defined as in Tab. [H 
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0^ r+r- ^ e + /i, 30 fb"^ 


Ma [GeV] 


200 


250 


Ns 


87.8 


136.7 


^cxp 


6.4 X 10-3 


1.1 X 10-2 




0.262 


0.412 


AM^/M^ [%] 


2.8 


3.5 



Table 4: Required number of signal events, Ns, with £ = 30 fb for a 5 a discovery in the 
channel (p ^ t^t^ —>■ e + fi. The other quantities are defined as in Tab. [H 

(see in particular Fig. 5.6 of Ref. [9] for the T~^r^ ^ /i + jet mode). 

Now we turn to the evaluation of the expected precision of the Higgs-boson mass mea- 
surement. In spite of the escaping neutrinos, the Higgs-boson mass can be reconstructed 
in the if, A — > rr channel from the visible r momenta (r jets) and the missing transverse 
energy, E^^^^, using the coUinearity approximation for neutrinos from highly boosted r's. In 
the investigated region of Ma and tan (3 the two states A and H are nearly mass-degenerate. 
For most values of the other MSSM parameters the mass difference of A and H is much 
smaller than the achievable mass resolution. In this case the difference in reconstructing the 
A or the H will have no relevant effect on the achievable accuracy in the mass determina- 
tion. In some regions of the MSSM parameter space, however, a sizable splitting between 
Ma and Mh can occur even for Ma ^ Mz- We will discuss below the prospects in scenarios 
where the splitting between Ma and Mh is relatively large. The precision AM^/M^ shown 
in Tabs. [1] - H] is derived for the border of the parameter space in which a 5 cr discovery 
can be claimed, i.e. with Ns observed Higgs events. The statistical accuracy of the mass 
measurement has been evaluated via 

AM^ ^ RaU 

A higher precision can be achieved if more than Ns events are observed. The corresponding 
estimate for the precision is obtained by replacing Ns in eq. ( |T71) by the number of observed 
signal events, N^^. It should be noted that the prospective accuracy obtained from eq. (ITTI) 
does not take into account the uncertainties of the jet and missing energy scales. In 
the r^r^ —>■ jets mode these effects can lead to an additional 3% uncertainty in the mass 
measurement [57]. A more dedicated procedure of the mass measurement from the signal 
plus background data still has to be developed in the experimental analysis. However, we do 
not expect that the additional uncertainties will considerably degrade the accuracy of the 
Higgs boson mass measurement as calculated with eq. f|T7j) . 
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4 Results 



The results quoted in Sect. [3] for the required number of signal events depend only on the 
Higgs-boson mass, i.e. the event kinematics, but are independent of any specific MSSM 
scenario. In order to determine the 5 a discovery contours in the M^-tan (3 plane these 
results have to be confronted with the MSSM predictions. The number of signal events, A'^ev, 
for a given parameter point is evaluated via 

N^^ = o-fefe^ X BR((/) r+r") x BR^^ x ^exp • (18) 

Here C denotes the luminosity collected with the CMS detector, a^^i^ is the Higgs-boson pro- 
duction cross section, BR(0 t^t^) is the branching ratio of the Higgs boson to r leptons, 
BRt-t- is the product of the branching ratios of the two r leptons into their respective final 
state, 

BR(r^jet + X) ^ 0.65, (19) 
BR(r ^/i + X) ^BR(r^e + X) ^ 0.175, (20) 

and £exp denotes the total experimental selection efficiency for the respective process (as 
given in Tabs. [H-Hl). The Higgs-boson production cross sections and decay branching ratios 
have been evaluated with FeynHiggs as described in Sect. 12. 2[ 

4.1 Discovery reach for heavy neutral MSSM Higgs bosons 

The number of signal events, A'ev, in the MSSM depends besides the parameters Ma and 
tan /3, which govern the MSSM Higgs sector at lowest order, in principle also on all other 
MSSM parameters. In the following we analyze how stable the results for the 5a discovery 
contours in the M^-tan (3 plane are with respect to variations of the other MSSM parameters. 
We take into account both effects from higher-order corrections, as discussed in Sect. 12. 2[ 
and from decays of the heavy Higgs bosons into supersymmetric particles. As starting point 
of our analysis we use the m™*^^ and no-mixing benchmark scenarios, where we investigate 
in detail the sensitivity of the discovery contours with respect to variations of the parameter 
fi. We then discuss the possible impact of varying other MSSM parameters. 

We have evaluated iVev in the two benchmark scenarios as a function of Ma and tan (5. 
For fixed Ma we have varied tan/3 such that N^^ = Ns (as given in Tabs. [H-H]). This tan/3 
value is then identified as the point on the 5 a discovery contour corresponding to the chosen 
value of Ma- In this way we have determined the 5 a discovery contours for the m™^^ and 
the no-mixing scenarios for /i = ±200, ±1000 GeV. 

In Figs. [T]-[3]we show the 5cr discovery contours obtained from the process btxp, t^t~ 
for the final states t~^t~ — jets, t~^t~ — > e + jet and r+r" jet. As can be seen 

from Tab. HI the fourth channel discussed above, r+r^ — > e ± /x, contributes for 30 fb~^ 
only in the region of relatively small Ma values and has a lower sensitivity than the other 
three channels. We therefore omit this channel in the following discussion. The discovery 
contours in Figs. [T]-[3]are given for the m™''^ and no- mixing benchmark scenarios with /i = 
±200, ±1000 GeV. As explained above, the 5 a discovery contours are affected by a change 
in fi in two ways. Higher-order contributions, in particular the ones associated with A;,, 
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Figure 1: Variation of the 5a discovery contours obtained from the channel bb(j), (j) — > t~^t — 
jets in the m^^^ (left) and no-mixing (right) benchmark scenarios for different values of /i. 
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Figure 2: Variation of the 5o" discovery contours obtained from the channel 660, (p ^ t~^t — >■ 
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Figure 3: Variation of the 5a discovery contours obtained from the channel bb(j), (p ^ t^t 
/i+ jet in the m^^^ (left) and no-mixing (right) benchmark scenarios for different values of /x. 



modify the Higgs-boson production cross sections and decay branching ratios. Furthermore 
the mass eigenvalues of the charginos and neutralinos vary with /i, possibly opening up the 
decay channels of the Higgs bosons to supersymmetric particles, which reduces the branching 
ratio to r leptons. 

The results for the 5 a discovery contours for the final state r+r^ jets are shown in 
Fig. [T]for the m™*^^ (left) and the no-mixing (right) scenario. As expected from the discussion 
of the Afe corrections in Sect. 12.21 the variation of the 5 a discovery contours with ix is more 
pronounced in the m™^^ scenario, where a shift up to Atan/5 = 12 can be observed for 
Ma = 800 GeV. For low Ma values (corresponding also to lower tan (3 values on the discovery 
contours) the variation stays below A tan /3 = 3. In the no-mixing scenario the variation does 
not exceed A tan/? = 5. The r+r~ jets channel has also been discussed in Ref. [14]. Our 
results, which are based on the latest CMS studies using full simulation [57], are qualitatively 
in good agreement with Ref. [14], in which the earlier CMS studies of Refs. [23,24] had beed 
used. The Sex discovery regions are largest for /i = —1000 GeV and pushed to highest tan/3 
values for ^ = +200 GeV. In the low Ma region our discovery contours are very similar 
to those obtained in Ref. [14]. In the high Ma region. Ma ~ 800 GeV, corresponding to 
larger values of tan (3 on the discovery contours, our improved evaluation of the 5 a discovery 
contours gives rise to a shift towards higher tan/3 values compared to Ref. [14] of about 
Atan/3 = 8 (mostly due to the up-to-date experimental input). Accordingly, we find a 
smaller discovery region compared to Ref. [14] and therefore an enlarged "LHC wedge" 
region where only the light CP-even MSSM Higgs boson can be detected at the 5 a level. 

The results for the channel r+r^ — > e + jet are shown in Fig. [2J Again the m™'^^ scenario 
shows a stronger variation than the no-mixing scenario. The resulting shift in tan /3 reaches 
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up to A tan /? = 8 for Ma = 500 GeV in the m^^^ scenario, but stays below A tan (3 = A 
for the no-mixing scenario. Finally in Fig. [3] the results for the channel r+r~ ^ + jet are 
depicted. The level of variation of the 5 a discovery contours is the same as for the e + jet 
final statelll 




Figure 4: Variation of the 5cr discovery contours obtained from the channel 660, r+r^ 
jets in the mf^^^ (left) and no-mixing (right) benchmark scenarios for different values of in 
the case where no decays of the heavy Higgs bosons into supersymmetric particles are taken 
into account (see text). 



In order to gain a better understanding of how sensitively the discovery contours in the 
M4-tan [3 plane depend on the chosen SUSY scenario, it is useful to separately investigate the 
different effects caused by varying the parameter /i. For simplicity, we restrict the following 
discussion to the 660, r^r^ jets channel. In Fig. H] we show the same results as 
in Fig. [H but for the case where no decays of the heavy Higgs bosons into supersymmetric 
particles are taken into account. As a consequence, the variation of the 5 a discovery contours 
with fi shown in Fig. H] is purely an effect of higher-order corrections, predominantly those 
entering via A;,. The difference between Fig. [T]and Fig. HI on the other hand, is purely an 
effect of the change in BR(0 r+r^) caused by the variation of the partial Higgs-boson 
decay widths into supersymmetric particles arising from a shift in the masses of the charginos 
and neutralinos. 

In Fig. m the dependence of the 5 a discovery contours on fi significantly differs from the 
case of Fig. [H While in Fig. [1] the inclusion of decays into supersymmetric particles gives 

^Since the results of the experimental simulation for this channel are available only for two Ma values, 
the interpolation is a straight line. This may result in a slightly larger uncertainty of the results shown in 
Fig. [3] compared to the other figures. 
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rise to the fact that the smallest discovery region is found for small /i values, /i = +200 GeV 
(with the exception of the region of very small Ma) , in Fig. H] the 5 a discovery contours 
are ordered monotonously in /x: the largest (smallest) 5 a discovery regions are obtained for 
/i = — (+)1000 GeV, i.e. for the largest (smallest) values of the bottom Yukawa coupling. 
As expected, the effect of the higher-order corrections is largest in the high tan/3-region 
(corresponding to large values of Ma on the discovery contours). In this region the variation 
of /i shifts the discovery contours by up to Atan/3 = 11 for the case of the m™'^'' scenario 
(left plot of Fig. H]), i.e. the effect is about the same as for the case where decays into 
supersymmetric particles are included. For lower values of tan (3 (corresponding to smaller 
values of Ma on the discovery contours), on the other hand, the modification of the Higgs 
branching ratio as a consequence of decays into supersymmetric particles yields the dominant 
effect on the 5 a discovery contours. Accordingly, the observed variation with n in this 
region is significantly smaller in Fig. H] as compared to the full result of Fig. [H The reduced 
sensitivity of the discovery contours on fi can also clearly be seen for the case of the no- 
mixing scenario (right plot), where as discussed above the A;, correction is smaller than in 
the m™^^ scenario. 




Figure 5: Variation of the 5cr discovery contours obtained from the channel bb(j), (f) ^ t^t^ 
jets in the m™'^^ (left) and no-mixing (right) benchmark scenarios with /i = +1000 GeV for 
different values of rrig. 

A parameter affecting the A;, corrections, see eq. (fTOj) . but not the kinematics of the 
Higgs-boson decays is the gluino mass, rrig. We now investigate the impact of varying this 
parameter, which is normally fixed to the values rrig = 800, 1600 GeV in the m™*^^ and 
no-mixing benchmark scenarios, respectively. The results for four different values of the 
gluino mass, nig = 200, 500, 1000, 2000 GeV, are shown in Fig. [51 The parameter has been 
set to yU = +1000 GeV in Fig. [5], such that the Higgs decay channels into charginos and 
neutralinos are suppressed. As one can see from eq. ffTOj) . the change of affects the 0{as) 
part of Afe and corresponds to a monotonous increase of A^. As an example, this yields for 
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/i = 1000 GeV, tan /5 = 50 in the two scenarios: 



max 



m~g = 200 GeV 



A, 



0.50 



h 



max 



m-g = 2000 GeV 
m-g = 200 GeV 
m-g = 2000 GeV 



A, 
A, 
Ab 



0.94 



no- mixing, 



0.06 



no-mixmg, 



0.29 . 



(21) 



In the no-mixing scenario the At value is close to zero, suppressing the mg-independent 
contribution to Af,, while the higher SUSY mass scale results in an overall reduction of Af, in 
this scenario. The value of A^ in the no-mixing scenario would slightly increase if rrig were 
raised to even larger values, but this effect would not change the qualitative behaviour. 

Fig. [5] shows that the results for the discovery reach in the M^-tan (3 plane are relatively 
stable with respect to variations of the gluino mass. The shift in the discovery contours 
remains below about Atan/3 = 4 for the m™'^^ scenario (left plot) and A tan/? = 1 for the 
no- mixing scenario (right plot). For the positive sign of fi chosen in Fig. [5], where the A^ 
correction yields a suppression of the bottom Yukawa coupling, the largest discovery reach 
is obtained for small rrig, while the smallest discovery reach is obtained for large rrig. This 
behaviour would be reversed by a change of sign of /i. 

We have also investigated the possible impact of other MSSM parameters (besides n and 
rrig) on the 5 a discovery contours in the M^-tan (3 plane. The A;, corrections depend also on 
the parameters in the stop and sbottom sector, see eq. ffTOj) . While the formulas in Sect. 12.2.21 
have been given for the region where MgusY ^ ^t, the qualitative effect of reducing the 
stop and sbottom masses can nevertheless be inferred. Sizable Ab corrections require relative 
large values of /i and rrig. If these parameters are kept large while the stop and sbottom 
masses are reduced, the A;, corrections tend to decrease. It is obvious from eq. (fTOj) that 
reducing the absolute value of At decreases the electroweak part of the A^ correction. As 
discussed above, this effect of the A^ corrections manifests itself in the comparison of the 
m™^^ and no-mixing scenarios, see Figs. [IHSl Concerning the possible impact of the Af, 
corrections on the 5 a discovery contours for the bb(j), — >■ r+r~ channel in the M^-tan f3 
plane we conclude that larger effects than those shown in Figs. [TH5] (where we have displayed 
the discovery contours up to tan /3 = 50) would only arise if the variation of fi were extended 
over an even wider interval than —1000 GeV < < -1-1000 GeV as done in our analysis 
above. 

We now turn to the possible effects of other higher-order corrections beyond those entering 
via Af, on the 5 a discovery contours for the bb(j), t^t~ channel. These effects are in 
general non- negligible, see the discussions in Sect. 12.21 and in Sect. 14.21 below, but smaller 
than those induced by A;,. As a consequence, the impact on the 5 a discovery contours in the 
Myi-tan (3 plane of other supersymmetric parameters entering via higher-order corrections is 
in general much smaller than the effect of varying /z in the high-tan (3 region of Fig. |H As 
an example, the difference observed in Figs. [TH5] between the rn^^^ and no- mixing scenarios 
arising from the different values of At and Msusy in the two scenarios (see eqs. f|T^ . f|T5l) ) 
is mainly an effect of the Af, corrections, while the impact of other higher-order corrections 
involving At and Msusy is found to be small. 

Also the decays of the heavy neutral MSSM Higgs bosons into supersymmetric particles 
are in general affected by other supersymmetric parameters in addition to the dependence 
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on /i, Ma and tan/?. The resulting effects on BR(0 t^t~) turn out to be rather small, 
however. We find that sizable deviations from the values of BR(0 —>■ t^t~) occurring in 
the m^^^ and no-mixing scenarios for —1000 GeV < < +1000 GeV are only possible in 
quite extreme regions of the MSSM parameter space that are already highly constrained by 
existing experimental data. 

Our discussion above has been given in the context of the MSSM with real parameters. 
Since the sensitivity of the 5 a discovery contours in the M^-tan /? plane on the other super- 
symmetric parameters can mainly be understood as an effect of higher-order corrections to 
the bottom Yukawa coupling and of the kinematics of Higgs-boson decays into supersymmet- 
ric particles, no qualitative changes of our results are expected for the case where complex 
phases are taken into account. 



4.2 Higgs-boson mass precision 

The discussion in the previous section shows that the prospective discovery reach of the 
bb(j), (J) — >■ T^T~ channel in the M^-tan j3 plane is rather stable with respect to variations of the 
other MSSM parameters. We now turn to the second part of our analysis and investigate the 
expected statistical precision of the Higgs-boson mass measurement. The expected statistical 
precision is evaluated as described in Sect. [3l see eq. ( |T7I) . In Figs. [6]-[7]we show the expected 
precision for the mass measurement achievable from the channel bb(j), t^t~ using the 
final states r+r~ jets and r+r" — >■ e + jet. Within the 5 a discovery region we have 
indicated contour lines corresponding to different values of the expected precision, AM/M. 
The results are shown in the m™*^^ benchmark scenario for = —200 GeV (left plots) and 
/i = +200 GeV (right plots). 

We find that experimental precisions of AM^/M^ of 1-4% are reachable within the dis- 
covery region. A better precision is reached for larger tan (3 and smaller Ma as a consequence 
of the higher number of signal events in this region. The other scenarios and other values of 
/i discussed above yield qualitatively similar results to those shown in Figs. [6], [3 

As discussed above, for large values of Ma the heavy neutral MSSM Higgs bosons are 
nearly mass-degenerate, Mh ~ Ma- The experimental separation of the two states H and A 
(or the corresponding mass eigenstates in the CP- violating case) will therefore be challenging. 
The results shown in Figs. |6] - [7] have been obtained using the combined sample of H and 
A events. It is important to note, however, that even in the region of large Ma the mass 
splitting between Mh and Mh can reach the level of a few %. An example of such a scenario 
is (as above, we consider the CP-conserving case, i.e. the MSSM with real parameters; the 
corresponding scenario in the case of non-vanishing complex phases has been discussed in 
Ref. [22]) 

MsusY = 500 GeV, At = A^ = 1000 GeV, /i = 1000 GeV, 

Ma = 500 GeV, Mi = 250 GeV, m~g = 500 GeV . (22) 

In Fig. [H] the mass splitting 

^Mha ^ \Mh-Ma\ . . 

M ~ mm{MH,MA) 

is given as a function of Xt for tan (3 = 40 and two Ma values. Ma = 300 GeV (solid 
line) and Ma = 500 GeV (dashed line). The dot-dashed and dotted parts of the contours for 
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Figure 6: The statistical precision of the Higgs-boson mass measurement achievable from 
the channel bbcf), cf) — > t+t~ — > jets in the m™^ benchmark scenario for fj, — —200 GeV (left) 
and fj, — +200 GeV (right) is shown together with the 5 a discovery contour. 
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Figure 7: The statistical precision of the Higgs-boson mass measurement achievable from 
the channel bb(f), (p t^t~ — e + jet in the m™*^^ benchmark scenario for /j, — —200 GeV 
(left) and /i = +200 GeV (right) is shown together with the 5 a discovery contour. 
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Figure 8: The mass splitting between the heavy neutral MSSM Higgs bosons, /S.Mha/M = 
\Mh — Ma \ /m.in[MH, M4), is shown as a function of Xt for Ma = 300, 500 GeV in a scenario 
with MsusY = 500 GeV, /i = 1000 GeV and tan/5 = 40. The other parameters are given in 
eq. ([22]). The dot-dashed (dotted) parts of the contours for Ma = 300 GeV {Ma = 500 GeV) 
indicate parameter combinations that are excluded by the search for the light CP-even Higgs 
boson of the MSSM at LEP [3]. 



Ma = 300, 500 GeV, respectively, in the region of small \Xt \ indicate parameter combinations 
that result in relatively low Mh values that are excluded by the search for the light CP-even 
Higgs boson of the MSSM at LEP [3]. One can see in Fig. [8] that the mass splitting between 
Mh and Ma shows a pronounced dependence on Xt in this scenario. Mass differences of up 
to 5% are possible for large Xt (while the widths of the Higgs bosons are at the 1-1.5% level 
in this parameter region). 

The example of Fig. [8] shows that a precise mass measurement at the LHC may in 
favourable regions of the MSSM parameter space open the exciting possibility to distin- 
guish between the signals of H and A production. In confronting Fig. [8] with the expected 
accuracies obtained in Figs. [6] - [7| one of course needs to take into account that a separate 
treatment of the H and A channels in Figs. [6|-[7|would reduce the number of signal events by 
a factor of 2, resulting in a degradation of the expected accuracies (for the same luminosity) 
by a factor of a/2. A more detailed analysis of the potential for experimentally resolving two 
mass peaks would furthermore have to include effects arising from overlapping Higgs signals. 
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Such an analysis goes beyond the scope of the present paper. 

5 Conclusions 

We have analyzed the reach of the CMS experiment with 30 or 60 fb^^ for the heavy neutral 
MSSM Higgs bosons, depending on tan/3 and the Higgs-boson mass scale, Ma- We have 
focused on the channel bbH/A, H/A — > t~^t~ with the r's subsequently decaying to jets 
and/or leptons. The experimental analysis, yielding the number of events needed for a 
5 a discovery (depending on the mass of the Higgs boson) was performed with full CMS 
detector simulation and reconstruction for the final states of di-r-lepton decays. The events 
were generated with PYTHIA. 

The experimental analysis has been combined with predictions for the Higgs-boson masses, 
production processes and decay channels obtained with the code FeynHiggs, taking into ac- 
count all relevant higher-order corrections as well as possible decays of the heavy Higgs 
bosons into supersymmetric particles. We have analyzed the sensitivity of the 5 a discov- 
ery contours in the M^-tan (3 plane to variations of the other supersymmetric parameters. 
We have shown that the discovery contours are relatively stable with respect to the im- 
pact of additional parameters. The biggest effects, resulting from higher-order corrections 
to the bottom Yukawa coupling and from the kinematics of Higgs decays into charginos 
and neutralinos, are caused by varying the absolute value and the sign of the higgsino mass 
parameter /i. The corresponding shift in the 5 a discovery contours amounts up to about 
Atan/3 = 10. The effects of other contributions to the relation between the bottom-quark 
mass and the bottom Yukawa coupling, arising from the gluino mass and the parameters in 
the stop and sbottom sector, are in general smaller than the shifts induced by a variation 
of jj,. The same holds for the impact of higher-order contributions beyond the corrections to 
the bottom Yukawa coupling and for the possible effects of other decay modes of the heavy 
Higgs bosons into supersymmetric particles. The results of our analysis, which was carried 
out in the framework of the CP-conserving MSSM, should not be substantially affected by 
the inclusion of complex phases of the soft-breaking parameters. 

We have analyzed the prospective accuracy of the mass measurement of the heavy neu- 
tral MSSM Higgs bosons in the channel bbH/A, H/A — > t'^t~. We find that statistical 
experimental precisions of 1-4% are reachable within the discovery region. These results, 
obtained from a simple estimate of the prospective accuracies, arc not expected to consid- 
erably degrade if further uncertainties related to background effects and jet and missing 
scales are taken into account. We have pointed out that a %-level precision of the mass 
measurements could in favourable regions of the MSSM parameter allow to experimentally 
resolve the signals of the two heavy MSSM Higgs bosons. 
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